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We propose a quasiprobability function, which is operationally defined and commensurate with 
experimental circumstance in the sense that its functionals exactly coincide with what they are 
supposed to represent. This enables to directly compare quantum statistics to those of a local hidden 
variable model with noninvasive measurability imposed on, called a classical model. We prove that 
the commensurate quasiprobability admits negative value(s) for a quantum state, whereas it is 
nonnegative for the classical model, being a legitimate probability function. It is remarkable that 
"nonclassicality" of a given quantum state is operationally determined. In addition we derive a 
sufficient condition for the entanglement of qudits in terms of a marginal quasiprobability of spatial 
correlations. 

INTRODUCTION 

Quantum physics contains distinct features from classical physics such as complementarity and nonlocality. It has 
generally been acknowledged how difficult to understand quantum physics and to unravel its features as they (or some 
of them) are counter-intuitive. As a transparent methodology, however, the comparison of quantum and classical 
statistics has provided significant intuitions in understanding quantum physics and separating its features from the 
classical. For instance, photons have been shown to occasionally exhibit antibunching that no classical statistics of 
light can imply. Those features are said "nonclassical" in a sense that the classical theory of light does not predict. 

To compare quantum with classical statistics, Wigner function has been employed on phase space, which corresponds 
to a probability distribution function in classical statistics Due to complementary (or uncertainty) principle, 

however, the Wigner function is not always positive semidefinite and negative- valued for some quantum states. As it is 
not allowed by any classical distribution, the negativity is regarded as a signature of the nonclassicality. The Wigner 
function, called a quasiprobability (distribution) function, has been generalized to discrete systems particularly of 
qubits as quantum informatics attains its importance [3J H] . The generalized quasiprobability has been applied to 
omnidirectional range of quantum information processing, including quantum tomography, quantum teleportation, 
and analysis of quantum algorithms [5j |6] . 

Even though they have established significant progress in their own context, the previous approaches are "in- 
commensurate" with experimental circumstances in the sense that all functionals do not necessarily coincide with 
what they seem to stand for. For instance, consider a functional of continuous- variable Wigner function W{x,p), 
C\W] = J dxdpxpW(x,p). One might regard this as the coincident or consecutive expectation, (xp) of position x 
and momentum p. This is de facto a different measurement whose basis is a set of the eigenvectors of a Hermitian 
observable operator (xp + px) /2 [2 ; those eigenvectors have nothing to do with the coincident nor consecutive mea- 
surements of the observable operators x and p. In the sense, any tries to directly compare such an incommensurate 
quasiprobability to its classical counterpart are pointless or misleading. In fact, the previous approaches make it 
difficult to interpret "classicality" [7J. 

In this Letter, we propose an alternative approach to define a commensurate quasiprobability function, by which 
we can directly compare quantum statistics to the classical. For this purpose it is crucial to specify a classical model 
and as one we take a local hidden variable model with noninvasive measurability imposed on [81110). We prove that 
for the classical model every commensurate quasiprobability is positive semidefinite, being a legitimate probability 
function. Based on the proof, we classify classical and nonclassical states. It is remarkable that the nonclassicality is 
operationally determined, depending on a given experimental circumstance. In addition we derive a sufficient condition 
for the entanglement of two qudits in terms of a marginal quasiprobability. 

COMMENSURATE QUASIPROBABILITY 
Quantum mechanical description 

Measurements are at the heart of quantum theory. Their description is also central in our defining a commensurate 
quasiprobability. Suppose that K possible (incompatible) observables are selectively and consecutively measured 
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for a given system. Each nondestructive and nondegenerate measurement for Ak is configured at a distinct space-time 
point, Xk = (tk, ffe) with t% < t% < ■ ■ ■ < tn, if it is to be measured. Every local measurement for Ak is assumed to 
be represented by a basis set {|a)fc} for D possible values, a £ {0, 1, D — 1}. For each run, we perform one of 2 K 
(consectutive) measurements, depending on the selection of the K observables, denoted by a tuple (tlx, n 2 , Uk ); 
Ak is selected if ^ 0. 

In quantum theory, we need to carefully describe consecutive measurements. For instance, suppose that two 
observables A\^ are selected to be measured. The local measurement for A\ collapses a given initial quantum state 
p to one of eigenstates |ai)i in the probability p(di) = i(o-i\p\ai)i- No local measurement is configured for A?.. That 
for A 3 collapses now the state |<Zi)i, previously collapsed into, to |a 3 ) 3 in the probability p(a 3 \ai) — | 3 (a 3 |ai) 1 1 2 . 
This completes the given consecutive measurement. We employ a specific form of expectations, 0, n 3 , 0, ■ • • ) = 

Sai~a3=o ui niai+n3a " 3 p{a 3 \a\)p{ai) , where ui — e 27Tl / D is a primitive D-th root of unity and n 13 are integers in the 
interval [1,D). All expectations of such a form compose a function, which we call a characteristic function and 
denote by x( n ) with n = (m, n2, uk). When local measurements are not restricted to the von Neumann projective 
measures, we can employ positive operator valued measures [llj . In the case each local measurement for Ak is 
represented by a set of Kraus operators {^(a)} with an overcompleteness relation J2 a =o ^l( a )^4fc( a ) = 1) where 1 
is an identity operator on the Hilbert space fid of dimension d. The characteristic function of a state p is given by 
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'[p], (1) 



(n) = TrTf] (1 - S ni , ) ^ iJ nia Ai(a) 



where T is an identity superoperator, T[p\ = p, and Ak(a) are Kraus superoperators for observable A^, Ak{o)[p\ — 
Aki^pA^a). The product of two superoperators is defined by their composition, for instance ^2-4i[p] = ^42^1/5^41^2' 
and T is a time ordering operator such that TAkAi = TAiAk = AkAi for tk > tj. Note that Eq. ([T]) results from all 
the possible (consecutive) measurements. 

We now define a commensurate quasiprobability function by a discrete Fourier transformation from x( n ) : 



™( a )^^f>~ a ' n x( n )< ( 2 ) 



n=0 



where a = (&i, a2, ax), the sum over n is the abbreviation of those over nk for k = 1, K, and the inner product 
a • n = Ylk=i a k n k- By definition, w(&) reproduces all the expectations x( n ) by the functionals, 
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X(n) = 5>™ W (a), (3 ) 

a=0 

where the sum over a stands for those over all ak- In other words, every functional of w(a) as in Eq. coincides 
with what it is supposed to represent, i.e., the expectation on the right hand side of Eq. |7p that results faithfully 
according to the experimental circumstance. The quasiprobability u>(a) is a real- valued function on the manifold of 
the K incompatible observables, which satisfies (Q.l) unit measure: ^] a ro(a) = 1, (Q.2) marginals: Summing w(sl) 
over a part of arguments results in a marginal quasiprobability of the rest, and (Q.3) the quantum probability for a 
single observable: The marginal quasiprobability of a single argument a^, Wk(ak) is the quantum probability of the 
appropriate measurement for Ak, that is, Wk(ak) = k{ a k\pWk)k- The importance of Q.2 and Q.3 has been emphasized 
for quantum tomography [HIS]. The commensurate quasiprobability can be negative- valued like Wigner functions, 
whereas its positivity is ensured by a certain class of classical theories. 



Local hidden variable theories with noninvasive measurability 



As a classical model, we take a local hidden variable model with noninvasive measurability throughout this Letter. 
Various types of hidden variable models have widely been adopted in their own context [8HTTH [T^HTS] . Nonlocality has 
been contrasted with local hidden variable models [8j [9l [13] . The local hidden variable models assume the existence 
of nonnegative joint probabilities involving all possible observations [9j[T3]. Our classical model also begins with such 
joint probabilities, assuming the noninvasive measurability that the values of observables are not altered by past nor 
future (whatever) observations |10) . This is understood as a not only spatially local but also temporally local hidden 
variable model, that we call the classical model. 
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FIG. 1. Degrees of nonclassicality N for pure states, corresponding to points on the Bloch sphere, for given observables (a) 



A% t 2 such that a X}V = (— l) Al < 2 and (b) Ai,2,3 such that a x ,y,z = ( — 1) 



The degrees are converted into the hue color 



system that the dark blue color implies the lack of the nonclassicality. The maximal degrees A/" m ax = (v2 — l)/4 is achieved if 
the qubit is in one of (a) four states p = (±l,±l,0)/v2 and (b) eight additional states p = (±1, 0, ±l)/v2, (0, ±1, ±l)/v2. 



The assumptions of our classical model replace the quantum characteristic function in Eq. ([I]) by the classical, 

D-l 

Xci(n) = £ c n > cl (a), (4) 



where p c \(a) are (nonnegative) joint probabilities on the outcomes of all the possible observables, determined by the 
hidden variables (see the discussions around Eq. (2) in Ref. [§])■ The quasiprobability function from Xci( n ) by the 
Fourier transformation, as in Eq. ([2]), is exactly equal to the probability function: u> c i(a) = Pel (a). This shows that 
in our classical model every quasiprobability is nonnegative and it is a legitimate probability function. 



Classical and nonclassical states 

Quantum theory conflicts with the classical model and it allows negative value(s) of a commensurate quasiproba- 
bility. If any, the negativity signifies "nonclassicality" in the sense that it can not be simulated by the classical model. 
A quantum state is said classical if its quasiprobability is nonnegative and otherwise nonclassical. 

When the local measurements are of mutually unbiased bases (MUBs) (in which D — d) [THUTT], the quasiprobability 
function in Eq. ^ becomes [T8] 



w(a) 




• p 



(•5) 
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where the generalized Bloch vectors p = TrAp and c?fe(a) = TrAAj,(a) with A = (Ai, A2, A^2_ 1 ). Here, {A^} is a 
complete orthogonal basis of Hermitian operators in Hilbert-Schmidt space of bound operators such that Ao = 1 and 
TrAiAj = dSij for i,j = 0,...,d 2 — 1. Every Bloch vector p obeys p ■ f > —1, for all Bloch vectors r, due to the 
positivity of quantum states [13 HH] ■ 

To illustrate nonclassical states, we consider a simple system of a qubit {d — 2) on which two observables are 
measured (K = 2). The two observable operators A i 2 are assumed such that a x = (— l) Al and & y = (— 1) A2 , 
where <J X:V ^ Z are Pauli spin operators. A quantum state is p — ^(1 + p ■ A), where the Bloch vector \p\ < 1 and 
A = (AijAajAji) = (a x ,a y , a z ). Then, the quasiprobability w{a\,a,2) = [1 + {—l) ai p x + (— l) a2 p y ]/4. For a state of 
p orthogonal to the x-y plane by the observable operators, w{a\,a2) is always positive semidefmite and the state 
is classical. On the other hand, a state of p x + p y > 1 is nonclassical as w(l, 1) = (1 — p x — p y )/A < 0. In 
particular, if a state is pure and its p is on the x-y plane, one of w(ax,a2) is always negative except p = ±x,±y. 
We quantify the nonclassicality by defining the degree J\f as the absolute sum of all the negative values of w(&): 
I\f = \ ^ a (|tu(a)| — w(&)). The nonclassicality degrees Af for pure states are presented in Fig.Fna), where the states 
correspond to points on the Bloch sphere and their degrees are converted into the hue color system with dark blue 
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FIG. 2. A proposal for the optical experiment to construct quasiprobability function of a polarization qubit, where D ai ,a 2 are 
photon detectors and PBS are polarizing beam splitters of the horizontal/vertical (H/V) and diagonal/anti-diagonal (D/A). 
Each PBS is placed on the path or placed away, depending on the (consecutive) measurement (ni,^)- 



indicating AT = 0. Here, the maximal degree of nonclassicality is A/" max = (V%— 1)/4 for four states p = (±1, ±1, 0)/ V2. 
In Fig. |l|b), the degrees are similarly presented when measuring three mutually complementary observables ^.1,2,3 
such that & XtViZ = (-l) Al > 2 - 3 , in which case K = 3 and w(ai,a 2 ,a 3 ) = [1 + {-l) ai p x + (-l) a2 p y + (-l) a3 p z }/8. Now 
all the pure states are nonclassical except p — ±i, ±y, ±z. The maximal nonclassicality is the same as the case of the 
two observables A\^. and located at eight additional states p = (±1, 0, ±1)/V2, (0, ±1, ±l)/-\/2. It is interesting that 
the classical states for the three mutually complementary observables ^1,2,3 coincide with those in Ref. l7j, while they 
do not in other cases. 

We would emphasize that the nonclassicality or classicality of a quantum state is operationally determined in our 
approach of the commensurate quasiprobability \17\ 119) . For instance, in the circumstance of measuring just two 
observables Ai 2 , certain states are classical whereas they becomes nonclassical in that of measuring Ai 2 3 (compare 
Figs, [lja) and '(b)). 



OPTICAL SYSTEM WITH INEFFICIENT DETECTIONS 



We propose an optical experiment of a polarization qubit to test if the classical or the quantum model properly pre- 
dicts an experimental result. The experiment of measuring two observables is schematically presented in Fig.[2j where 
D ai ,a 2 are photon detectors and PBS is a polarizing beam splitter of the horizontal/vertical (H/V) or diagonal/anti- 
diagonal (D/A). Each PBS is placed on the path or placed away, depending on the (consecutive) measurement (ni, n 2 ) 
of two observables, where = 0, 1. 

To obtain an experimental quasiprobability w exp , letting Ai j2 be observables of H/V and D/A, we construct each 
expectation of (ni,n 2 ) from the set of counts {N niri2 (ai, 02)} at the detectors D ail a 2 and the total count number 
N ni n 2 - Xoxp(«i,«2) = J2a u a 2 ujniai+n2a2 f ni n 2 {ai,a 2 ), where / ni „ 2 (ai,a 2 ) = N niTl2 (at, a 2 )/N niri2 . Then, w cxp (a 1 ,a 2 ) 
is given by the Fourier transformation as in Eq. (pi). 



SUFFICIENT CONDITION FOR ENTANGLEMENT 



When applied to a composite system consisting of spatially separated subsystems, its commensurate quasiprobability 
reveals two types of nonclassicality. One is of the temporal correlation(s) of each subsystem which we discussed 
previously. The other is of the spatial correlations of the subsystems, that originates from entanglement. The 
segregation of the latter from the former can be made by considering a marginal quasiprobability of spatial correlations, 
which is nonnegative for any separable state whereas negative for an entangled state. 

To illustrate, we consider two qudits 1 and 2 in a state p on the Hilbert space Hd <S> Hd with d a prime or power of 
prime. Each qudit j is measured on K{= d + 1) local observables Aj k , as for the single qudit. The quasiprobability 
u>(a) is a function of a = (ajk), where ajk are the outcomes of Ajk- We take w c (c) = 53a(Ilfc ^c k .ai k +a 2k ) w(a) as a 
marginal quasiprobability of spatial correlations. Here the Kronecker delta S a jj = 1 if a and b are congruent modulo 
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d. The marginal quasiprobability is given |18j by 

Wc (c) = ^ T (l + TrSM(c)), (6) 

where S = Tr A (8 X T p and \ T = (\J , AT 2 _,) with the transposition T defined in the standard basis {!«)}• Here, 

M ( c ) = ^y^y^"2,fc(cfc - a)ai,fc(a), (7) 

k a 

where the generalized Bloch vectors cfife(a) = TrXAik(a) and c?2fc(o) = TrA T j4.2fc(a) for the two complete sets of MUB 
{Ajk(a)} of qudits j = 1, 2. For each c, the matrix M(c) maps a generalized Bloch vector of qudit 1 to that of qudit 
2, preserving the vector norm |18) . 

For an arbitrary separable state p = YliPiPu ® (hi, S = Y^iPiPufki and Eq. (44) becomes 



Note that p^i ■ M(c) ■ pu > — 1 for each i, because M • pu is a Bloch vector of qudit 2 and every pair of Bloch vectors 
have an overlap no less than —1 due to the state positivity [17j[T8]. Thus, the marginal quasiprobability is positive 
semidefinite for every separable state, whereas it can be negative for some entangled states jTSj. This reciprocally 
implies the sufficient condition that a quantum state is entangled if its marginal quasiprobability of spatial correlations 
has any negative value. In the sense, the marginal play a role of entanglement witness |20j . 



SUMMARY 



we proposed a commensurate quasiprobability function to faithfully describe (consecutive) measurements in the 
operational way. It was proved that it is positive semidefinite and a legitimate probability function for every local 
hidden variable model with the noninvasive measurability. We explicitly illustrated classical and nonclassical states for 
a qubit and showed that the nonclassicality is operationally determined, depending on give experimental circumstance. 
In analyzing the entanglement we proposed a marginal quasiprobability of spatial correlations to segregate spatial 
nonclassicality from temporal nonclassicality and showed that its negative value, if any, is a sufficient condition for 
two-qudit entanglement. It will be interesting to apply to quantum information processing, for instance, testing if 
any found algorithm is classically simulated, and to Bell's inequalities for an arbitrary system. It is an open question 
how to define a commensurate quasiprobability for a continuous variable system, which might suffer from theoretical 
difficulties when its observables are unbound [21] . 
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SUPPLEMENTARY INFORMATION 



We derive explicitly the crucial results presented in the paper. For the purpose we employ the generalized Bloch 
representation and provide some properties in the notations of the paper. For a single qudit, we derive the form of 
commensurate quasiprobability function when local measurements are given by mutually unbiased bases. For two 
qudits, we derive a marginal quasiprobability of spatial correlations, in addition to the characteristic and quasiprob- 
ability functions. We prove the sufficient condition for two-qudit entanglement that a quantum state is entangled if 
its marginal quasiprobability of spatial correlations has any negative value. As a nontrivial example of applying the 
theorem, we classify one-parameter family of two-qudit Werner states. 




Generalized Bloch representation 

A generalized Bloch representation is convenient and employed to derive the crucial results in the paper. We 
summarize the results given by Ref. [17] and provide some additional properties in the notations of the paper. 

Consider an orthogonal complete set of Hermitian operators, {Ai\i = 0, d 2 — 1}, on the Hilbert-Schmidt space of 
bound operators such that Ao = 1 and Tr A^Aj = d5ij where 5ij is a Kronecker delta. By using the set, a Hcrmtian 
operator H and a density operator p are expanded to 

d 2 -i \ 

(9) 
(10) 

where hi = TrA^-ff and pj = Tr Xip. For every density operator (unit-trace positive operator) p, po = 1 and we call 
p= (pi, Pd 2 -i) a generalized Bloch vector (GBV). We employ a concise notation hereafter: 

p = Tr A/3, (11) 

where A = (Ai, A2, AcP-i)- Because Trp 2 < 1, the norm of p is upper bounded: \p\ < \/d — 1. If p is pure, 
\p\ = V d — 1 . The GBVs stay within a Bloch sphere of radius \/d — 1. However, not all vectors within the Bloch 
sphere correspond to density operators, that is, there is no one-to-one correspondence between density operators and 
vectors within the Bloch sphere. In fact, the set of GBVs specifying density operators is restricted by the positivity 
of density operators such that a given density operator p should hold 

Trpr > p-f > -1, (12) 

for all pure density operators f with |r| = \fd — 1. 

A density operator for two qudits, p is given in the Bloch representation by 



1 - - d 1 - ~ \ 

p=-p\t(g>l + p 1 -\(g)l + t(g>P2-\ T + ^ S <*P>>a ® J , (13) 

where pj are local GBVs of qudit j and S = Tr A ® X T p. Here A T = (X[ , Aj , Aj 2 _ x ) and the operator transposition 

A„ is defined in the standard basis {|n)}. We call S a generalized Bloch matrix (GBM). From the upper bound 

s2 



Tip 2 < 1, the local GBVs and GBM obey, 



The positivity of p is represented by 



J2 \Pj\ 2 +TrS T S < d 2 - 1. (14) 

3=1,2 



}2 P,-r J +TrST T >-l, (15) 

3=1,2 
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where fj are the local GBVs and T is the GBM of an arbitrary pure state f . If the local GBVs p\ = p^ = such as 
maximally entangled states and Werner states, the inequality ( 14 ) becomes 

TrS T S<d 2 -l (16) 



and the inequality (15 1 of state positivity becomes 

TrST T >-l. (17) 

In particular, if we are focused on any pure product state f = fi (£> Tjj, then T = f\fi so that the state positivity 
implies 

fi-S-f a >-l. (18) 
Some properties of GBVs and GBMs are summarized. 
(G.l) A vector p is a GBV if and only if satisfying the two conditions; 



(D.l) being within the Bloch sphere, \p\ < \J d — 1, and 

(D.2) the state positivity as in Eq. (12), p ' ■ f > —1 for every pure GBV f with |r| = \/d — 1. 

(G.2) Given a GBV p, the vector a = £p is a GBV for < f < 1: 

It is clear that a satisfies the condition (D.l). It also does (D.2) as 

ff.f = £p-r > -£> -1, (19) 

for every pure GBV f. 

(G.3) A convex combination of GBVs is a GBV: 

Let pi be GBVs and < ^ < 1 with = 1. The vector p = ^ £iPi satisfies (D.l) as its norm is no larger 

than the largest among the GBVs pi and thus \p\ 2 < d — 1. The vector p also satisfies (D.2): For every GBV f, 



p ■ T 



= -t>-]T6 = -1, (20) 



where we used pi ■ f > — 1 . 
(D.3) A function from the set of GBVs to itself is called a Bloch map. 

(G.4) For a two-qudit state p with local GBVs pj = 0, its GBM S = Tr A (8) X T p satisfies f\ ■ S > —1 for an arbitrary 
pair of GBVs fj, as shown in Eq. (18 1. 

(G.5) An orthogonal GBM S of a two-qudit state p is a Bloch map from qudit 2 to qudit 1: 

For each GBV f, the image a = S ■ f satisfies the condition (D.l) as S is orthogonal. As TrS T S = d 2 — 1, 
by Eq. (14) all local GBVs of p vanish. By (G.4), the image a satisfies (D.2). Thus S is a Bloch map by the 
definition (D.3). 

A nondegenerate and orthogonal measurement A is represented by a set of basis projectors {A(a) = \a)(a\}. The 
basis projectors are unit-trace positive operators, satisfying the orthonormality relation 

8 a , al = TrA(a)A{a') = ^(1 + 3(a) ■ S{a% (21) 



where a (a) are GBVs of A(a), a(a) = Tr XA(a). Eq. (21 ) is rewritten as a(a) ■ a(a') = dS a ^ a ' — 1. The basis projectors 
also satisfy the completeness relation, given as 



d-l 



^A(o) = !<*=>■ 5(a) = 0. (22) 

a=0 a=0 

Their GBVs for the measurement A, a(a) form an overcompleteness relation in the (d — l)-dimensional subspace of 
Bloch sphere: 



^ d-l 

- ^ a(a)a(a) = ld-i, (23) 



a=0 



9 



where a(a)a(af) is a tensor product and is an identity matrix on the (d — 1) -dimensional subspace. 

For a qubit with d — 2, a nondcgcneratc and orthogonal measurement A is given by two projectors, say A{a) — \a) (a\ 
for a = 0, 1. In the Bloch representation they are given as 

A{a) = l(l + (-l) a a-X), (24) 



where a is a unit vector of "orientation" for A and A = (a x , a y ,a z ) is the vector of Pauli spin operators, a x ,y,z- The 
GBVs of projectors A(a) are thus a(a) — (— l) a a. 

Two measurements and A 2 are mutually complementary if precise knowledge in one of them implies that all 
possible outcomes in the other are equally probable J22j [23] , 



Tr A 1 (a)A 2 (b) = -, Va, b = 0,...,<2- 1 (25) 

where {Ak(a) = |a) fcfc (a| } is the set of basis projectors for A^. This leads to the orthogonality relations of 

ai(ai) ■ d 2 (a2) = 0, Vai,a 2 = 0, (26) 

where <3fc(a) are GBVs of the basis projectors A^{a) for A^. It is then notable that the two mutually complementary 
measurements (MCM) have two exclusive Bloch subspaces, each of (d— 1) dimension. As the Bloch space is (d 2 — 1)- 
dimensional, there exist at most (d+ 1) MCM. The existence of (d+ 1) MCM was proved for prime or power of prime 
dimensions |16j . Even though the proof does not exclude the possibility of other dimensions, we assume that our 
subsystems are prime or power-of-prime dimensional to avoid any unnecessary confusion. The set of all MCM is said 
complete. The complete set of MCM allows the overcompleteness relation in the Bloch sphere, 

^ d+1 d-1 d+1 

2 ^k{a)a k (a) = ^2 lfe,d_i = ld 2 -i, (27) 

fc=la=0 k=l 

where l^-.d-i is an identity matrix on the fc-th exclusive Bloch subspace and Icp-i is that on the Bloch sphere. The 
basis sets of MCM are said mutually unbiased bases (MUB). 
We add another property of GBVs related with MUB. 

(G.6) For the given set {<Sfc(a)}, the two operator sets of {A Q } and {Ak(a)} are related by 

d+1 d-1 

X = Y / J2 d k(a)Ak(a) (28) 

fc=l a=0 

A k (a) = 3 (l + a k {a)-\) . (29) 



The first transformation results from Eqs. (22| and (27), and the second does from Eqs. (26), (21), and (22). 

(G.7) Quantum tomographic equivalence: 

Consider a complete set of MUB, {Ak(a)} with the set of GBVs {cSfe(a) — Tr XAk(a)}. The complete set expands 
a quantum state p as 



P = 



, / d+1 d-1 \ 

-(a + ^^ft(fl)4(«)] (30) 

k=la=0 



d+1 d-1 

^5> fc (a)i fe (a)-l, (31) 

k=l a=0 



wher e w e used the transformation as in Eq. (28) and pk(a) = p ■ dk(a) = dpk(a) — 1 with Pfe(a) = Tr Ak(a)p. 
Eq. (31) shows that every quantum state is completely determined by the complete set of MUB {Ak(a)} and 
their probability distributions pk{a) (or the coefficients pk(a)). Consider another state, defined as 

, / d+1 d-1 \ 

\ fe=la=0 / 

d+1 d-1 

= 5m^( a )^( a )- a - ( 33 ) 



k=l a=0 
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where the set {pk(a) = dpk{a) — 1} is the same as that of p and {Bk(a)} is another complete set of MUB. We 
say that a is equivalent to p in the sense that both are in the same form with respect to their complete sets of 
MUB and in particular both have the common probability distributions pk(a) — Bk{a)a = TrAk(a)p. The 
two operators a and p are connected by the relation, 

^ d+l d-l 

3 = Pk(a)a k {a)-p, (34) 

fc=l a=0 

where /3fc(a) = TrA-Bfe(a). This implies that the matrix M = h ^ fe ^ a /3k(a)ctk(a) is a Bloch map from one 
complete set of MUB {a>k(a)} to another {/^(a)} (see (D.3)). As an example, consider a state p which is 



diagonal in the basis {Ai(a)}: p — ^2 a P\(.a)Ai(a). This is equal to Eq. (31 1 with pk(a) = 1/d for all pairs of a 



and k 1. The operator a = ^2 a Pi(a.)Bi(a) is a quantum state and a = ^ ^ Q /3i(a)(? 1 (a) • p is a GBV. 



Quasiprobability for a qudit 



We shall derive the form of commensurate quasiprobability function w(sl) when local measurements are of mutually 
unbiased bases (MUB). The quasiprobability function in Eq. (2) of the paper is rewritten as 



w (a) = Tr 



K 



. k=l 



D 



AA(a/c) [p] 



(35) 



where a = (oi, 02, ctif), A-/4fc(&) = Ak{a) — j> X) ' Ak{a') for each fc, and p is a quantum state. Here I is an 
identity superoperator, I[p] = p and Ak(a) are Kraus superoperators of measurement Ak, Ak(a)[p] — Ak(a)pA' k (a). 
The product of two superoperators is defined by their composition, for instance _42-4i[/5] = A2AipA\At 2 , and T is a 
(time) ordering operator such that TAkAi = TAiAk = AkAi for k > I. 

If Ak are of MUB, they satisfy mutual complementarity in Eq. ( 25 ), each with D = d outcomes, and so two different 
measurements Ak and Ai satisfy 



TtA k (a)Ai(a')[p] = \ r TtAi{a!) 

d 



(36) 



for each a and a'. This leads to Tr AAk(a)AAi(a')[p] = 0, V/5. The generalization is straightforward to arbitrary 
combinations of different measurements of MUB: For each (ak, a;), 



Tr AA k (a k ) ■ ■ ■ AAi(ai)[p] = 0, V/5. 
In the case, the quasiprobability function is given as 



w(a) = 



^Tr (l + dJ2 ^4fc(ofc)J [p] = ]g f 1 + E « fc ( afc ) ' ^ > 



(37) 



(38) 



where p is the GBV of state p and ctk(a) are those of projectors Ak(a). Eq. (38) is the result presented in the paper. 
For a qubit, Eq. (38 1 is reduced to 



f(a)= ^ ^l+X;(-l)°*5 fc -^ , 



(39) 



where o?fc is the unit vector of the orientation Ak (see Eq. (24)). 



Quasiprobability of two qudits 



Consider a composite system consisting of two spatially separated qudits in a quantum state p. Each qudit is to be 
selectively and consecutively measured by a complete set of K = d + 1 MUB with D — d, as described in the paper 
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for a single qudit. The two complete sets of MUB are not necessarily the same. Its characteristic and quasiprobability 
functions are respectively given, similarly to Eqs. (1) and (2) in the paper, as 



2 d+l 

X (n) = Tr(g)rn 

j=l k=l L 



d-l 



S njk , l + (1 - S njkt0 ) J2 ^ ka A jk (a) 



a=0 



and 



2 d+l , v 

a) = Tr(g)T JJ I - + AA jk (a jk ) J [p], 

A — 1 I 1 V / 



(40) 



(41) 



j=i k=i 



where a = (o,-fc) with its elements being the value ajk of fc-th observable Ajk for qudit j and AAjk(a) — Ajk(a) — 
2 Sa' ■Ajk(a') for each (j, fc) (the other symbols are defined similarly to Eqs. (1) and (2) in the paper). By the similar 
reasoning to Eqs. p7| and (38), the quasiprobability results in 



^ / d+l d+l d+l \ 

d 2(d+l) I 1 + ^2 ®lk{aik) ■ "2i( a 2i) ■ 92 + ^2 "lfe( a lfe) • S • "2i( a 2i) , 

' y k=l 1=1 k,l=l J 



(42) 



where the local GBVs pi = TrA (g> Ip and p 2 = Tr 1 (g) A T p, the GBM S = Tr A <8> A T p, ai, fe (a) = Tr AAi jfe (a), and 
"2,fc(a) = Tr A T i 2 ,fe(a)- 



Marginal quasiprobability of spatial correlations and two-qudit entanglement 

We define a marginal quasiprobability of spatial correlations by 



d-l /d+l \ 

=( c ) = n <Wifc+«« ) w ( & )' 

a=0 \/c = l / 



(43) 



where c = (ci, c 2 , Q2+1) and the Kronecker delta £ ai f, = 1 if a and b are congruent modulo d and <5 a! b = otherwise. 
Here the sum over a is the abbreviation of those over ajk for j = 1, 2 and fc = 1, 2, d+l. By the relation for djk(a) 
as in Eq. (22), the marginal becomes 



w c (c) 

where the trace is defined on the Bloch sphere and 



d d - 



-(l + TrSM(c)), 



d+l d-l 



M(c) = ^ Y X! "2,fe(cfe - a)ai,fc(a). 



(44) 



(45) 



fe=l a=0 



Here, we note that the argument Ck — a is defined by its positive residue modulo d. 

The matrix M(c) for each c is a Bloch map from qudit 2 to qudit 1, in other words, it maps a GBV of qudit 2 to 
qudit 1 (see (D.3), (G.5), and (G.7)). The matrix M(c) is an isometry (or orthogonal), preserving the vector norm: 
Letting M = M(c) for a given c, 



-1 d-l 



M T M = — ^ X! ^,k{a)d 2 ^k{ck ~ a) ■ a 2 ,;(q - 6)dij(6) 



k,l=l a,b=0 
d+l d-l 



\^2^2ai.A a )aiA c 



k = l a=0 
ll,d 2 -l, 



(46) 



where we used the relations as in Eqs. (21), (22), (26), and (27). Here, l^-i is an identity matrix on the Bloch 
space of qudit j. Similarly, 



^ d+l d-l 

MM T = - ^ "2,fc(cfc - a)a 2 ,k{ck - a) = $-2,<P-i- 



(47) 



k=l o=0 
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These properties of M(c) are employed in the paper to prove the sufficient condition that a quantum state is entangled 
if its marginal quasiprobability of spatial correlations w c (c) has any negative value. 

The marginal function iw c (c) in Eq. (44) plays a role of entanglement witness [20J . As a nontrivial example, we 
consider a Werner state of two qudits, p w — f\ip)(ip\ + (1 — f)t <g> t/d 2 , where \tp) = ^2 n \n,n)/Vd is a maximally 
entangled state and < / < 1. In the Bloch representation, it becomes 



d 2 -i 

1 ® 1 + / ^2 S al3 X a 

aP=l 




(48) 



where S = (-0|A <E) A T |-0) = 5 a p is the GBM of \ip). Note that the GBM S is orthogonal and it is a Bloch map from 
qudit 2 to qudit 1 so that for f>2 a GBV of qudit 2, p\ = S ■ p2 is a GBV of qudit 1 (see (G.5)). The marginal functon 
w c (c) is given from Eq. (44 1 as 



d'l- 



-(l + /TrSM(c)). 



This implies 



Tr SM(c) = ^J2Y1 dlk ^ ■ S ' a ^ Ck ^ a )^^EE = ^ +1 ) 



(49) 



(50) 



where we used (D.2). The marginal function is thus lower bounded, 

w c (c)>J^[l-(d+l)f] 



(51) 



The equality in Eq. (51 1 holds for c = (0, 0, ...0) if the complete set of MUB for qudit 2 is chosen such that a.2k{o) = 
S T • dik(ck — a + l) for each (k, a). Then, the minimum is negative as far as / > l/(d + 1), implying that the Werner 
state p w in Eq. (48 ) is entangled if / > l/(d + 1). 



